tt a tt St Ht tt HH HH HH Ha 
tt a tt a a tt i 


T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 


HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HH HHH 
HHHHHHHHHHHHHHH 
HHHHHHHHHHHHHHH 
HHHHHHHHHHHHHHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 


tt tt tt st 4 Ht Ht Ht Ht Ht Ht HH HH 
AAA a 
| LE) 


LLLLLLLLLLLLLLL 
LLLLLLLLLLLLLLL 
LLLLLLLLLLLLLLL 


**F JLE®*1D**MTHGSINH 


MM MM 
MM MM 
MMMM MMMM 
MMMM MMMM 
“4 6M UMM 
“Mm 66MM MM 
MM MM 
MM MM 
MM MM 
MM MM 
MM MM 
MM MM 
MM MM 
MM MM 
LL 
ae 
LL 
LL 
LL 
LL 
LL 
LL 
LL 
LL 
LL 
LL 
LLLLLLLLLL 


LLLLLLLLLL 


ast at tt st tt ot 
me a ee a ee eg ee eg 


— 
—e 
— ee 
Oo bees 


es Gem Sheet Geert Pema Gad fener feet) (wo Gene ed oe be 
Pa bad et ee ee Det ed Gt et ee ee 


oo 
| inne 
ee 
cio 


HH HH 
HH HH 
HH HH 
HH HH 
HH HH 
HH HH 
HHHHHHHHHH 
HHHHHHHHHH 
HH HH 
HH HH 
HH HH 
HH HH 
HH HH 
HH HH 
SSSSSSSS 
SSSSSSSS 
SS) 
SS 
SS 
$$ 
SSSSSS 
SSSSSS 
SS 
SS 
SS 
SS 
SSSSSSSS 
SSSSSSSS 


GGGGGG6GG 
GGGGGGGG 


GG GGGGGG 
GG lat 


GGGGGG 
GGGGGG 


eee 


os 
oe 


— 


— 
fe es et Cs 


be ne bs Bs me ms a be Pe a an ts 


_— 
ee 
ee 
os 


MTHS$GSINH 


Table of contents 


: 


i 


6 
; G Floating Hyperbolic Sine venlen 


HISTORY ; Detailed Current Edit Wistory 
DECLARATIONS  ; Declarative Part of Module 
MTHSGSINH = Standard G Floating DSINH 


16-SEP-1984 01:31:30 VAX/VMS Macro V04-00 


6 
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1-006 a aes BE MTHRTL.SRCIJMTHGSINH.MAR; 1 ° (1) 
-TITLE MTHSGSINH 3 F Cleasing Hyperbolic Sine routine 
. IDENT /1-006/ > File: MTHGSINH.MAR EDIT: JCW1006 


s Seneenueenceceeseeeeeeeneennonnenennseseesenenceeeoreneeeesceceeeerqeonecees 
. 

COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 

DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 

ALL RIGHTS RESERVED. 
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INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER * 
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OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
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DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
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3 LARRRRRAAASAAAASALARASASASZARERRRRERRRR RE RARER RRR R RR RRR ERR RRR R RRR RR ER ER RRR ERE 
> FACILITY: MATH LIBRARY 
be 
; ABSTRACT: 

MTHS$SGSINH is a function which returns the G floating yyeerbet te sine 


of its G floating point argument. The call is standar 
call-by-reference. 
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So 
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; VERSION: 1 

; HISTORY: 

; AUTHOR: 

Steven B. Lionel, 26-Jan-79: Version 1 


; MODIFIED BY: 
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; Detailed Current Edit History 6=-SEP-19 :24:16 CMTHRTL. SRC CM THGSINH.MAR; 1 1-0 
-SBTTL HISTORY ; Detailed Current Edit History 


ALGORITHMIC DIFFERENCES FROM FP-11/C ROUTINE: none 
Edit History for Version 1 of MTHS$GSINH 


WO OONOUS wr" 


SOOCOCOCOCOCOOCOOOOOCOOoOoOo S- 


SOOGCOCOCOCOCCOSCOCOOSOOoCOoOSO <~Oo 


$0 3 1- Gdeotee from + ha? hg serene, 1-001. SBL 26-Jan-79 

3 é : Ie 8 - Use MTHSGEXP ept 97 

33 ; ; 1- - eUaingted Sand ett to EXP a input values between 27*\n2 

: an 

38 61 ; = Changed all final floating point divisions by 2 to interger 
62 ; subrtracts of 1 from the exponent field. 

33 65 ; - Extended maximum range 1024*\n2. 

38 rt: ; - sneneges, he for computing EXP(ix!=ln2) to reduce error. | 

00 66 ; 1-004 = Changed W* to G* in call to MTHS$SIGNAL RNH gg ecaee -1981 

00 67 ; 1-005 = Eliminated eyes te short Literals. RNH 15 

00 68 ; 1-006 - Changed a MNEGD to a MNEGG. JCW 20-Jan-83 


MTHSGSINH 
1-006 


bet 


00000004 


So 
So 
So 
So 


20€ 28723 87084052 
7B606E3A 288740A6 
S9EEFEFA 2E4240A6 
SCOOFEFA 26424006 
4(67361C 86548050 


oouUnmmo—for-oOm 
ee ee nae 

VI — Pr OMOC 
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- SBTTL 


INCLUDE FILES: 
EXTERNAL SYMBOLS: 
DSA 


6 
erbolic Sine re. Sh 16-SEP-19 
eclarative Part of Modul 6-SE 


DECLARATIONS 


7 
3 F 
5g 
° ye 
7 : L GBL 
0 8 SEXTRN MTHSSSIGNAL 
28 SEXTRN MTHSK_FLOOVEMAT 
$ 44 -EXTRN MTHSGEXP_R6 
00 f + EQUATED SYMBOLS: 
00 5 value =4 
8 
fy § : MACROS: none 
000 87 ; 
000 88 
$00 «90 
909 31 > PSECT DECLARATIONS: 
000 98 .PSECT _MTHSCODE 
a8 
Mita 96 : OWN STORAGE: none 
8 
000 _99 : CONSTANTS: 
B33 190 
000 108 G_27_L0G_2: 
000 10 -QU “%20E 2872387084052 
008 104 G_1023_L0G 2: 
008 105 , “X7B606E 3A28B740A6 
01 106 6_1024_L0G_2: 
01 10 .Q0AD “X39EEFEFAZE4240A6 
01 108 G_LOG_2_HI: 
01 109 .QUAD “X3COOFEFAZE424006 
0 8 110 G_LOG_2_L0: 
0 1] -QUAD *X4C67361C86548D50 
6 11§ GSINHTAB: 
0 115 . WORD 20037172.20167 15 
¢ 116 -WORD “00 $146, 20147 0 
1 ~WORD 2003 347,°001674 
4 118 -wORD *00 544 +2012625 
8 11 “WORD *003 31 »*0 0640 
¢ 120 .WwORD *0 13 50,°0044164 
4 121 <WORD *0037641,°0 1042) 
+ Ve ; <wORD * 1942120 Q 5 
48 1 .wORD * 3 7 2-20 : 39 
4c 126 ~WORD  *0052525,°0052526 
09 125 .WORD 0. 


P=19 


4 01 
4 11 


31:99 


; value.rg.r 


PIC,SHR,LO™S ,EXE,NOWRT 
3 program section for math routines 


3 27*Ln2 
3: 1023*ln2 
3 1024*Ln2 


; (high 43 bits of Ln2)+2e2#-43 


; bn2 = G_LOG_2_HI1 


; DECIMAL: 


; DECIMAL: 0.00 


AX/VMS Macro V04-00 
4:16 (CMTHRTL.SRCIMTHGSINH.MAR; 1 


; Declarative Part of Module 


0.2508223608819151D-07 
; DECIMAL: 0.2755729803646086D-05 
> DECIMAL: 0.1984126984813681D-03 
; DECIMAL: 0.83333333333352327D-02 
; DECIMAL: 0.16666666666666670+00 


MTHS 


F 6 
it a beCanation® S’Bectarative Bart of Modul 'S-SEP=1986 1:24:16 EATMMTL.SRESRTMGSINN.maR:1 "29% (3), 110 


127 
00000006 O038 135 GSINHLEN = .> GSINHTAB/8 
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MTHSGS INH G Floating Hyperbolic Sine rontes 16-SEP-1984 01:31:30 VAX/VMS Macro v04-00 ® 5 TH 
1-006 ATHSGSINN 2 Standard G Flosting DSINH  GoSEP-}oBe S]iausim EMTANTe emesaTNogINNmaRs? Pe! dy 1-01 
] -SBTTL MTHSGSINH = Standard G Floating DSINH 
: 4 ;++ 
2 ! 5 ; FUNCTIONAL DESCRIPTION: 
05 5 ; GSINH = G floating point function 
B2 133 : GSINH(X) is computed as: 
058 141 ; If {Xi < 2ee-27) GSINH(X) = X. | 
3 148 3 IF Qe*-27 =< {Xi <_0,.25, GSINH(X) = Chebyshev Series. 
1435 ; If 5 =< ix} 2788s Ms baa = (GEXP(X) = GEXP(-x))/2, 
05 144 ; If 7* n2 =< {Ki < 1025*ln2, GSINH(X) = AT pad pet EEE D 
05 145 ; If 198 ps =< Mt < 1024*ln2, then GSINH(K) = sign(X)*GEXPCixXi-Ln2). 
93 146 ; If 1024*in2 =< 1X}, then overflow. 
058 148 ; CALLING SEQUENCE : | 
0058 130 ; GSINH.wg.v = MTHSGSINH(x.rg.r) 
0098 138 > INPUT PARAMETERS: | 
00000004 005 154 LONG = ; define Longuerd multiplier 
00000004 Bpe 133 x = 1 * LONG ; Contents of x is the argument | 
Ba 38 139 3 IMPLICIT INPUTS: none 
0098 159 > OUTPUT PARAMETERS: 
Ones 11 ; VALUE: G floating hyperbolic sine of the argument | 
0038 168 : IMPLICIT OUTPUTS: none | 
0058 165 ; COMPLETION CODES: none | 
0058 167 : SIDE EFFECTS: 
005 169 : Signal: MTH$_FLOOVEMAT if 1024*Ln2 < {Xi with reserved operand in RO/R1 | 
005 170 ; (copied to the signal mechanism vector CHFS$L_MCH_RO/R1 by LIBSSIGNAL). 
005 171 ; Associated message is: ‘FLOATING OVERFLOW IN MATA LIBRARY'’. Result is reserved 
bne 76 ; operand -0.0 unless a user supplied (or any) error handler changes CHFSL_MCH_RO/R1 | 
05 174 : NOTE: This procedure disab! floating point underflow, enables integer | 
05 175 ; overflow. 
05 176 : | 
05 177 3e-- 
gh 
40FC 005 180 ENTRY MTHSGSINH, “M<IV, R2, R3, R4, RS, Rb, R7> 
OSA 181 ; standard call-by-reference entry 
OSA 1 g ; disable DV (and FU), enable IV 
Ben 1 MTHSFLAG_JACKET ; flag that this is a jacket procedure in 
6D Q0000000'GF 9 03h MOVAB G*MTHSSJACKET_HND, (FP) 
61 ; set handler address to jacket 
0061 ; handler 
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MTHSGS INH ; G Floating Hyperbolic Sine routine 16-SEP-1964 01:71:30 VAX/VMS M v04-00 Pp MTH' 
1-006 ATHSGSINN 2 Standard G Flosting DSINH  GeSEP-19Be Siizasie EMTARTE.caedarncsiamaRs? 2% (fy 1-0 
1 
$1 184 ; case of an error in routine 
61 185 ; If an error, convert signal to user PC 
61 1 § : and resigna 
ce 9% c 50FD 6} MOVG ayalu (AP), R6 : R6/R7 = 1X! = avalue(AP) 
0 g000 FAA 069 1 § BICw2 #*x8000, R pA Ute IM! 
FFO BF «=650.—sé@I 06 190 CMPW = sCéRO, -#*XSFF t compare |X! with 0.25 
1€ 18 191 BGEQ  GEd_T0_0.2 ; branch if ix! >= 0. 
a 
75 194: IK! < 0.25 | 
Be 
3660 BF 50 B61 0075 19 CMPW 8 =s_-RO,-: #*X3E60 ; ‘Xt with 2-27 
04 18 OO7A 198 BGEQ GEG. T0.2M27 + gonnere il aT so Seed? | 
o7¢ 199 
o7C 00 ; 
O7C 201 : IK! < 200-27 
bore 508° | 
50 56 7D O07C 204 MOVQ = R6, _-RO : RO/R1 = X 
04 O07F 05 RET 3 return with result = argument 
$0 | 
0080 08 : 2ee-27 =< IK! < 0.25 
0080 540 | 
50 50 44FD a089 eit ee ee Get ARG**2 for POLYG | 
. bd ee j 
9f AF OOS 80 SEED 08s 518 POLYG RO. #GSINHLEN-1, GSINHTAB ne | 
OBA 214 ; RO/R1 = SUM(Ci*X**#i), with 
OBA 15 ; last coefficient zero 
50 56 44FD OOBA 216 MULG2 R6,RO ; MULG2 by ARG, and then 
50 56 40FD 008 17 ADDG R6,R0 ; add in ARG with overhang. 
04 009 18 RET : return with result in RO/R1 
nat 
09 1: 0.25 =< ix! 
093 958° 
009 4 GEQ_T0_0.25: 
FFOF CF 50 S1FD 009 5 CMPG RO, G_1023 96 2 3 compare {Xi with 1023*\n2 
3c 14 90 $ BGTR  GTR_TAAN_1023_[0G_2 > branch if {Xi > 1023*Ln2 
98 3 
tee : 3 0.25 =< IKt =< 1023*Ln2 
8 53h | 
FFSF CF 50 SIFD $98 § CMPG RO, G_27_L0G_2 ; Compare {Xi to 27*ln2. If greater 
20 = o14 in 7. BGTR ONLY_ONE- TERA z; only one call to GEXP is necessary. 
A 5 : 0.25 =< IK! < 27#Ln2 
AS 539° 
50 56 7D OOA 8 MOVQ = R6,_-RO ; RO/RI = xX 
00000000'EF 16 O0A6 JSB8 MTHSGEXP_R6 > RO/R1 = GEXP(X) 


6 
MTHSGS INH ; G Floating Hyperbolic Sine veuties 16-SEP-1984 01:31:30 \VAX/VMS Macro V04-00 Page 7. 
1-006 THSGS INH : Stendard G Floating DSINH ~SEP-1984 94:31:30 YRTHRTL. SREINTNGSINH.MAR: 1 ° (4). 
7E 50 7D OOAC 240 mova RO, =(SP) ; push GEXP(X) on stack 
38 086 BC S52FD AF 41 MNEGG @value(AP), RO : e /R1 = =x 
9 0 org 16 4 4¢ J MTHSGEXP_RO + RO/R1 = GEXP(-x) 
3 8 bag A 4 SUBG3 RO, (SP)#, RO ; RO/R1 = GEXP(X) = GEXP(-x) 
10 A F 264 SuBW = #*K0010, RO + RO/R1 = (GEXP(X)=GEXP(=x))/2 
04 ; c2 RET 3 return with result in RO/R1 
cS 57 
C 48 ; 27*ln2 =< IK! < 1023*Ln2 
c 559° | 
o¢ i ONLY_ONE_ TERM: 
00000000" EF 16 C ¢ J$B MTHSGEXP_R6 ; RO/R1 = GEXPC(IX!) 
04 OE ae OOD 58 pee. poeiitue™ ayy Ae 
3 negative change sign o 
50 50 52FD OOCF 555 MNEGG RO, RO S GEXPCTK!) ronkigis. | 
Op8 2$ POSITIVE: 
50 10 a Be 2h SUBW #*x0010, RO 3 RO/R1 = sign(X)*GEXP(iX!)/2 
thd 59 ; 
0007 60 ; 1023*ln2 =< {x} 
0007 61 3 
007 6¢ 
0D7 63 GTR_THAN 102s _LOG_2 
FF33 CF 50 SIFD 0007 264 TMP RO. G .1024_L0G6_2 : Compare 1X} to 1024*Ln2 
22. 18 Bpoe 22 Bees ERROR” 3: if 1024*ln2 =<ix!, overflow occurs 
OO 567 | 
OODF 68 ; 
BppE +4 : 1023" ln2 < IX! < 1024*Ln2 
ODF 0; 
OODF 71 
50 FF34 CF 42FD ODF Le: SUBG2 G_LOG 2 pil RO 3 RO/R1=!Ki-G_LOG 2 HI 
00000000'EF 16 O0E5 27 JSB MTHSGEXP Ré : RO/R1 = GEXP(!xXT-G_LOG_2_HI) 
52 50 FF30 CF 45FD oeB 74 MULG3 LOG 2 "CO: RO, R2 ; RO/RS=GENP(iX}~G_LO6_2“HI)*6_LOG_2_LO 
50 52 sere Fe 75 SUBG2 4 RO ~ 3 RO/RI=GEXPCiXi-Ln2) 
04 Bf 55FD OOF6 76 TSTG @value(AP) ; test the sign of X 
4 18 OOFA 77 BGEQ 10$ ; branch if X >= 0 
50 50 52FD th 78 MNEGG RO, RO ; RO/RI = si n(X) ® GEXPCIXi-bn2) 
04 0100 79 10$: RET 3 return with result in RO/R1 
101 280 
191 81 
191 28¢ | 
6104 
101 4 ; 1024*ln2 =< iXi, error 
191 S$; | 
7E OO°8F 9A 181 8S ERROR: nourel. #MTHSK_FLOOVEMAT, -(SP) ; condition value 
50 01 OF 79 0105 Hy #15, #T, RO 3 RO/ aT = result = reserved operand -0.0 
109 8 : goes © signal mechanism vector | 
109 90 3 chest OncH RO/R1) so error handler 
109 91 ; con modify the result. 
00000000'GF 01 #=FB 0109 3 CALLS #1, G*MTHSSSIGNAL : tgnet error and use real user's PC 
1 8 9 3 . ependent of CALL vs JSB 
04 BF 3 RET ; return - RO/R1 pestered from CHFSL_MCH_RO/ 
0111 96 


es Onderd Sine routlon,, Bee 8s P1iSkitg YATARTE. SaeSarncgiam.mars1 P29 dy 


Symbol table 


ERROR 101 R 1 
GEO. 10.0459 R 1 
GEQ_ TO 2m R 1 
GSIRHLEN = 

GSINHTAB R 1 
GTR_THAN. 1023_L0G_2 D7 R 1 
6.1 ¢7-L06. 08 R 1 
G_1024 LOG. 10 R 1 
G27 LOG 00 R 1 
et38-5- Sat 
LONG ~~ = 90000004 
MTHSSJACKET_HND teereere x 1 
MTHS$$S1GNAL eeeeeeee =X 0 
MTHSGEXP R6 eeeeeeee =X 0 
MTHSGSINA 00000058 RG i 
MTHSK_FLOOVEMAT eeeeeeee xX 0 
ONLY [ONE _TERM 000000C3 R 01 
POSITIVE 00000003 R 01 
VALUE = 00000004 
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4:16 CMTHRTL. JMTHGSINH.MAR; 1 


! Psect synopsis ! 
Feeoeooeoecoocone + 


PSECT name Allocation 

» OB «= 00000000 <( q-? 00 ( 0. 

_MTHSCODE 00000111 ¢ 273.) O1¢ 1. 
+ 

Phase Page faults CPU Time 

Initialization 3 90:00:00 .07 

Command processing 10 0:00:00. § 

Pass 1 86 0:00:00. 

Symbol table sort 9 0:00:00.01 

Pass 2 6 :00:00.76 

Symbol table output 4 8:08: 8-83 

Psect synopsis output 0:00:00.01 

Cross-reference output 8:82: 9-28 

Assembler run totals 29 :00:02.4 


The working set Limit was 900 pages. 


1 page of virtual memory was used to define 1 macro. 


LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
RD NOWRT NOVEC LONG 


Oooooooooo 
Sooooooeoco 
—-DOoCOooo°oo 
—COOOrvoOwso 
se 8 ff #& &@ *& &© @ 

SoolLwoviw 
MONMNCO— fm 


4544 bytes (9 pages) of virtual memory were used to buffer the intermediate code 
There were 10 pages of symbol table space allocated to hold 


d 1 local symbols. 
source Lines were read in Pass 1, producing 11 object records in Pass 2 


21 non-local an 


er 8epa lobe V1i3ciie EMTARTe.eaeimrucsina.mar;1 29° 19, 
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; G Floating Hyperbolic Sine ro cine 
winett ne Macro Run Statistics ~ ’ 


teen ecerenne See eee noon wae awa + 


! Macro Library statistics ! 


tere ree ewe ee cw ene ewe eee were $ 


Macro Library name Macros defined 

“$255$DUA28: CSYSLIBISTARLET.MLB;2 0 

0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 

MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHGSINH/OBJ=OBJ$:MTHGSINH MSRC$:MTHJACKET/UPDATE=(ENH$:MTHJACKET) #MSRC 


MTH! 
Tab! 


L EQUIPMENT CORPORA 
ENTIAL. AND PROPRIET 
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VAX/VMS V4.0 
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